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2 : ABSTRACT 

■ The Universidad Complutense de Madrid (UCM) Survey is a long-term project aim- 
ing at finding and analysing star-forming galaxies detected by their Ha emission in 
Schmidt objective-prism plates. The instrumental set-up limits the volume of the Uni- 

^ ' verse surveyed to a redshift z ^ 0.045. So far we have discovered several hundred 

■ emission-line galaxies (ELGs) undergoing star formation at a wide range of levels. An- 
alyzing a complete sample of ELGs from the UCM survey we have computed the Ha 
luminosity function for the star-forming galaxies in the surveyed volume of the Universe. 
A Schechter function provides a good fit to the Ha luminosity function with the follow- 
ing parameters: a = -1.3 ±0.2; L*(Ha) = iQ^s-isio.oSg^gg-i ^* ^ 10-3-2±0-2 Mpc'^ 
for Hq = 50kms~^ Mpc~^. Integrating over the full range of luminosities, we obtain an 
Ha luminosity density of 10^^'^^°'^ ergs~^ Mpc~^. Using the Ha emission as a star for- 
mation rate (SFR) estimator, this translates into a SFR density for the Local Universe 
of O.O13tE];|j[}^M0yr-^Mpc-^ in star-forming galaxies with £^M^(Ha+[NH]) > lOA and 
z <^ 0.045, for a Scalo Initial Mass Function. This is the first observational determina- 
tion of this quantity, which will provide a direct test of current galaxy formation and 
evolution models. 

Subject headings: Galaxies: luminosity function, fundamental parameters 
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1. Introduction 

Studying the properties of star forming galaxies at 
low redshift, and quantifying the amount and loca- 
tion of star formation in the Local Universe is a nec- 
essary step towards understanding galaxy formation 
and evolution. It can provide a crucial test for galaxy 
formation and evolution models. Such models are now 
able to put together cosmology, dark matter, gas and 
stars, and make definite predictions for the properties 
of the present-day galaxy population and its evolution 
with redshift. In particular, they can predict current 
SFRs (Kauffmann, White & Guiderdoni 1993, Cole 
et al. 1994) but very few observational constrains ex- 
ist. Substantial data sets have been put together (e.g. 
Kennicutt 1992, and references therein) but not for a 
star-formation selected sample of galaxies. 

Moreover, extensive observational effort has been 
devoted to the study of the properties and evolution of 
high redshift galaxies (e.g. deep counts, deep redshift 
surveys, gravitational lensing) but one of the major 
problems that arises when analysing high-z data is 
how to make a meaningful comparison with local sam- 
ples, usually obtained using very different observa- 
tional techniques. Carrying out such a comparison in- 
volves to distinct steps: first we need well-understood 
data for local samples that are directly related to the 
distant ones, and second, we have to transform the ap- 
parent properties of the local objects to predict how 
they would look like at high- 2; (i.e., introduce cosmo- 
logical corrections and selection effects). If cither or 
both of these steps is not accurate enough, apparent 
differences between local and distant samples could 
mimic evohition, and lead to erroneous conclusions. 
A clear example where this problem is very important 
occurs in the interpretation of the so-called 'faint blue 
galaxies' problem, where much of the recent contro- 
versy depends very much on what different authors 
believe about the properties of local galaxy samples, 
and in particular about the local galaxy LF. See, e.g., 
Koo & Kron (1992) and Broadhurst, Ellis & Glaze- 
brook (1992) for different views on this problem. 

Closely related to this is the recent evidence found 
in deep galaxy redshift surveys selected in the op- 
tical and in the near infrared which suggests star- 
formation activity substantially increases with red- 
shift from z ~ to 2; ~ 1 (Songaila et al. 1994, 
CoUess 1995, Ellis et al. 1995, Lilly et al. 1995). Al- 
though the causes of this increase in star formation 
are not yet clear (interactions and mergers have been 



suggested), it could explain the faint blue galaxy ex- 
cess and the suggested evolution in the optical field 
LF (see Aragon-Salamanca 1995 for a recent review). 
Paradoxically, we scxnn to know more about the num- 
bers of star-forming galaxies at intermediate and high 
redshifts than locally. This letter tackles this problem 
directly, providing the first published estimate of the 
numbers of nearby star-forming galaxies and the star 
formation rate (SFR) density of the Local Universe. 

The Ha luminosity is a very good direct measure- 
ment of the current star formation rate, since it is 
directly related to the number of massive stars (see, 
e.g., Kennicutt 1992, Gallagher and Gibson 1993). 
It is better than other optical Balmer lines such as 
H/3 or H7, affected more by stellar absorption and 
reddening. Metallic nebular lines like [OIII] A5007, 
[Oil] A3727 (affected by excitation and metallicity) 
or IRAS fluxes (affected by the dust abundance and 
properties) are rather star-formation indicators than 
quantitative measurements (see Gallagher et al. 1989, 
Kennicutt, 1992). Although it could be argued that 
near infrared recombination lines such as Br7 could be 
even better measurements of the current SFR, being 
less affected by extinction than the optical lines, the 
lack of near-IR instrumentation with very large field 
of view makes large area near-IR surveys impractical 
at the moment. These considerations imply that the 
best available way to quantify current star formation 
in the local Universe is by using an Ha-selected sam- 
ple of galaxies. However, it should be kept in mind 
that if some fraction of newly formed stars is com- 
pletely hidden by dust, they would not contribute to 
the Ha luminosity, implying that the SFR derived 
here should be considered a lower limit. 

The Universidad Complutense de Madrid (UCM) 
Survey for emission-line objects provides an ideal tool 
for such purpose, since galaxies were selected ac- 
cording to their Ha emission. It is a low-dispersion 
objective-prism survey for low-redshift emission-line 
galaxies (ELGs) being carried out using the 80/120 cm 
f/3 Schmidt telescope at the German-Spanish Obser- 
vatory of Calar Alto ( Almeri'a, Spain) . A 4° full aper- 
ture prism, yielding a dispersion of 1950 A/mm, com- 
bined with IIIaF emulsion, has been used to search 
for ELGs selected by the presence of Ha emission in 
their spectra. The instrumental setup is able to record 
the Ha-|-[NII] AA6548,6584 blend in emission for ob- 
jects up to z <^ 0.045. The survey itself and the two 
first discovery lists have been described in detail by 
Zamorano et al. 1994 & 1995. Preliminary results 
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and details are presented in Rego et al. (1989, 1993), 
and Zamorano et al. (1990, 1992). So far, the UCM 
Survey has found 264 candidates in 17 fields cover- 
ing 471.4 square degrees. The overall object density 
is 0.56 candidates per square degree. More than half 
of the sample (138 candidates) are galaxies which do 
not appear in any published catalogue. Follow-up 
imaging and spectroscopy of the survey candidates 
has been completed and used to confirm the candi- 
dates as ELGs and classify them. Morphologically, 
the UCM sample represents a heterogeneous popula- 
tion of galaxies dominated by late-type galaxies (66% 
being Sb or later) with about 10% presenting typical 
parameters of earlier types, and another ~ 10% be- 
ing irregulars. In all but 22% of the galaxies we find 
the emission to come almost entirely from the nuclear 
regions. Their median luminosity is about one mag- 
nitude fainter than that of magnitude-limited galaxy 
samples selected in blue plates, with a higher fraction 
of low-luminosity galaxies. As classified from their 
spectra and IRAS data, the most commonly found 
ELGs (47%) are intermediate to low-luminosity ob- 
jects with a very intense star-formation region which 
dominates the energy output of the galaxy. Their 
metallicities range from solar values to -^Zq, peak- 
ing at jZq. a full analysis has been published in 
Vitores et al. (1995a, 1995b & 1995c) and Gallego 
et al. (1995a, 1995b & 1995c). This extensive data 
set provides the raw materials for the study presented 
here. 

From this catalogue, we have built a representa- 
tive complete sample of star-forming galaxies (which 
excludes AGNs) suitable for the statistical study pre- 
sented in this Letter. In section 2 we discuss the 
galaxy selection and estimate the Ha luminosity func- 
tion, which in section 3 is used to estimate the present 
day SFR density. Finally, we summarise our conclu- 
sions in section 4. 

2. The Ha luminosity function 

Direct information on the amount and distribution 
of the present-day SFR in the Local Universe can 
be obtained by constructing the luminosity function 
(LF) for galaxies with current star-formation activ- 
ity. The Balmer line emission, and in particular the 
Ha luminosity, provide a direct signature of the star- 
formation activity occurring in the galaxies. We will 
now estimate the Ha LF. 

In a sample built from an objective-prism survey 



such as the UCM one, the survey completeness de- 
pends primarily on the line flux, not the apparent 
magnitude. As shown by Salzer (1989), the com- 
pleteness of a line-selected sample depends on the 
emission-line-f continuum flux and the contrast be- 
tween the line and continuum (the equivalent width 
of the line). Applying the procedures of Salzer (1989) 
to the UCM sample (Gallego 1995c) we deflne an ar- 
bitrary magnitude scale ttil-i-c as a function of the 
linc+continuum flux -FL+c (in ergs^^cm^^): 

mL+c - -17.0 - 2.5 log(F^L+c) (1) 

A quantitative estimate of the survey completeness 
can be obtained using the V/Vmax test (Schmidt 1968, 
Huchra and Sargent 1973). Such a test indicates that 
all the galaxies in the UCM sample with mL+c < 17.3 
and EW(Ra+[mi]) > lOA have U/Kiax ^ 0.5, thus 
constituting a complete sample (see Gallego et al. 
1995c for details). For the remainder of this study, the 
sample of 176 UCM galaxies that fulfil these criteria 
will be referred to as "the UCM representative com- 
plete sample". The ttil+c < 17.3 hmit corresponds to 
a line plus contimmm flux of 1.9 x 10^^'* erg s~^ cm~^. 

With such a sample we are able to compute the 
Ha luminosity function. The total Ha luminosity, 
reddening-corrected using the Balmer decrement, has 
been taken from Gallego et al. (1995b). Distances 
have been estimated from the galaxy redshifts, as- 
suming a value of Ho — 50kms^^ Mpc^^. Note that 
when considering the object detection procedure, the 
combined Ha4-[NII] emission has been taken into ac- 
count, since the objective-prism produces very low 
resolution. However, when computing Ha luminosi- 
ties, the [Nil] contribution is not included since the 
follow-up spectroscopy separates both lines. 

If $[logL(Ha)] is the number of galaxies per unit 
volume per 0.4 interval in log i(Ha) , then the Schmidt 
estimator can be written as 

$[logL(Ha)] = ^^-^ (2) 

where D, is the surveyed solid angle and V^^x represent 
the volume enclosed by a sphere of radius equal to the 
maximum distance the galaxy could have and still 
be detected in the survey. The summation is over 
all galaxies in the Ha luminosity range log L (Ha) ± 
0.5AlogL(Ha). We have used AlogL(Ha) = 0.4. 
The luminosity function for the UCM galaxies is given 
in Tabic |l|. 
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The Vrnax method presented here has the advan- 
tage of simphcity, but because galaxies are assumed 
to be distributed homogeneously, the results may be 
distorted at the faint end if there are local inhomo- 
geneities in the sample. The maximum likelihood esti- 
mator described in Mobasher, Sharpies & Ellis (1993) 
is independent of such inhomogeneities, provided that 
the LF has the same shape everywhere. Using this 
method, we obtain a LF that agrees, within the er- 
rors, with the one obtained from the Vmax method, 
showing that local inhomogeneities are either non- 
important or average out in the volume sampled by 
the UCM survey. 

Typical Ha luminosities range from 10^° to 10^^ erg s 
The density of star-forming galaxies decreases with 
Ha luminosity, as pointed out previously by other 
authors (Salzer 1989, Boroson et al. 1993). How- 
ever, and what is more important, the Ha technique 
used for the UCM survey recovers a higher number of 
star-forming galaxies at lower relative levels of star- 
formation activity than surveys carried out in blue 
region of the spectrum. This selection effect has clear 
repercussion when estimating the number and distri- 
bution of star-forming galaxies from such blue sur- 
veys, since they can miss up to ~ 50% of the galaxies 
undergoing star formation and bias the Ha luminos- 
ity distribution (thus the SFR distribution) towards 
higher luminosities (SFRs). 

3. The SFR density of the Local Universe 

One of the most important applications of the Ha 
luminosity function is that it can be used to estimate 
the current "SFR function" of galaxies, which de- 
scribes the number of star-forming galaxies as a func- 
tion of their on-going SFR. Moreover, if an analytical 
expression is obtained, it can be used to extrapolate 
observed luminosity density to total luminosity den- 
sity in the volume considered. The luminosity func- 
tion for our complete sample of current star-forming 
galaxies is well fitted by a Schechter (1976) function 

(p{L) dL = (j)* {L/L*)'^ exp(-L/L*) d{L/L*) (3) 

where the luminosity function (j){L) is directly related 
to $[logi(Ha)] for L = Z/(Ha) by the equation 

^{logL)^^ = <f>{L)dL (4) 

The best-fitting parameters (excluding the first point) 
are 

a = -1.3 ±0.2 



^* = 10-3-2±n-2Mpc-3 

L* = 10^2.15±0.08g^gg-l 

Figure 1 shows the binned luminosity distribution, 
and the solid curve represents the Schechter function 
with the above parameters. 

The errors associated with the derived parameters 
were obtained using a Monte-Carlo method. A large 
number of simulations were computed with errors in 
the LF following a Gaussian distribution of a equal 
to the square root of the number of galaxies in each 
log L (Ha) bin. A Schechter function was then fitted 
to the simulated LFs, yielding a normal distribution 
for the parameter values. The quoted errors are la 
for these distributions. 

If a Schechter function with the above parame- 
ters is a good approximation to the luminosity dis- 
tribution in the volume of Universe considered, since 
a < —2, 4){L) is finite and can be integrated for whole 
range of luminosities. 

POO 

Ltot = / (t>{L)LdL 
Jo 

= (j)* L* r(2 + a) (5) 

The gamma function r(2 + a) takes the value 1.30 
for a = 1.3, yielding a total Ha luminosity per 
unit volume of lO"^^'^^"'^ ergs"^ Mpc^'^ in the Local 
Universe {z ^ 0.045) for star-forming galaxies with 
EW(Ha+[mi]) > lOA. 

We can now relate this Ha luminosity to the SFR 
following the method introduced by Kennicutt (1983). 
We use the evolutionary models of Chariot & Bruzual 
(1993) to relate the current SFR to the number of ion- 
ising photons produced by the newly formed stars, 
and then Case B recombination theory to predict 
the luminosity of the Ha emission line. A Scalo 
(1986) Initial Mass Function, including stars in the 
O.IM0 < M < 125Af0 mass range was used. The 
transformation can be written as 

L(Ra) = 9.40 x 10^°-^^^^^ erg s'^ (6) 

The top axis of Figure 1 shows how the L(Ha) 
scale transforms into a SFR scale using this equa- 
tion, so that the current "SFR function" can be 
directly read from the figure. The total Ha lu- 
minosity density translates into a SFR density of 
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0.013lo:oo5 Mq, Mpc^^ where the main source of 
uncertainty originates in the normahzation factor </>* 
of the LF. 

Some caution is necessary when interpreting these 
quantitative results. The Ha luminosity is only sen- 
sitive to the SFR in stars with M > 10 Mq, the main 
contributors to the ionising flux. Thus, the total 
SFR given here is an extrapolation assuming a rea- 
sonable (local neighbourhood) IMF. To transform the 
total SFR values given here into SFRs in stars with 
M > IOMq, a factor of 0.044 should be appUed. Nev- 
ertheless, we have decided to quote total SFRs since 
they can provide useful constraints for galaxy forma- 
tion and evolution models, notwithstanding the IMF 
uncertainty. 

4. Summary and Conclusions 

Using an Ha selected sample of nearby galaxies we 
have estimated the Ha luminosity function for the lo- 
cal Universe. We argue that, since the Ha emission 
provides a good estimate of the on-going star forma- 
tion, the galaxies have been selected by their current 
SFR. Therefore, we can use the Ha luminosity func- 
tion to determine the "SFR function" describing the 
number of star-forming galaxies as a function of their 
SFR. Integrating over all Ha luminosities (or SFRs) 
we determine the current SFR density of the local 
Universe for galaxies with iiJM^(Ha-|- [Nil]) > 10 A and 
z ^ 0.045. The transformation of Ha luminosities 
into SFRs depends on the stellar initial mass func- 
tion. The values presented in Section 3 correspond to 
a solar neighbourhood Scalo (1986) IMF, but can be 
easily obtained for any IMF by changing equation 6. 

We stress that this has been done for the fist time 
for a sample of star-forming galaxies directly selected, 
as far as it is possible today, by their SFR. Using 
Ha as the SFR estimator has clear advantages over 
optical surveys carried out in the blue and far infrared 
surveys. It provides a cleaner sample, less biased by 
metallicity, excitation and dust-properties effects, and 
covers a very broad range in star-formation activity, 
reaching relatively low levels. 

Both the SFR function and the total SFR den- 
sity determined here can provide direct constraints 
to galaxy formation and evolution models. They can 
also play an important role in studies of high redshift 
galaxy samples providing the z = benchmark to be 
compared with the numbers of high-^; star-forming 
objects. 
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Table 1 

Ha Luminosity Function for the UCM Survey 



log h{Ha) 
[erg s-i] 


log(4>) 
[# per Mpc'"' per 0.4 
iiitcrN'Ml of logfLflln)] 


Number of 
galaxies 


40.6 


-3.380 


7 


41.0 


-2.878 


26 


41.4 


-3.081 


43 


41.8 


-3.323 


54 


42.2 


-3.654 


36 


42.6 


-4.564 


10 
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Fig. 1. Ha luminosity density distribntion for a 
complete sample of UCM galaxies. The errors are the 
square root of the number of galaxies in each bin. The 
solid line represents the fitted Schechter function to 
the data points. The scale at the top of the diagram 
indicates the SFR corresponding to the Ha luminosity 
(see equation 6). 
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